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In this brief report, we show the equivalence between the tight-binding descriptions of the mono- 
layer and bilayer honeycomb lattices. With appropriate value of the third nearest neighbors cou- 
pling, the Hamiltonian for a monolayer is equivalent to the low energy effective Hamiltonian for 
bilayer in the presence of trigonal warping. A simple physical argument is provided to explain this 
correspondance. 
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INTRODUCTION 



II. MONOLAYER WITH THIRD NEIGHBOR 
COUPLING 



Although the literature on physical properties of mono- 
layer and bilayer graphene based on the tight-binding 
description of the electronic spectrum is quite large, 1 it 
may be useful to precise a few simple relations between 
these two systems. In this brief report, we consider the 
tight-binding Hamiltonian on the monolayer honeycomb 
lattice with coupling between nearest neighbors (t) and 
third nearest neighbors (£3). It was noticed that for a 
critical value of the coupling ts (= t/2), the low en- 
ergy spectrum is quadratic, similar to the case of bilayer 
graphene. 2 It turns out that the low energy spectrum 
in bilayer s is not exactly quadratic, due to a small addi- 
tional hopping contribution between the layers leading to 
a trigonal warping of the spectrum. 3 Here we show that 
close to the critical value of the coupling £3 (< t/2), there 
is an exact correspondance with the low energy Hamil- 
tonian for bilayer graphene in the presence of trigonal 
warping. We show that there is also a simple correspon- 
dance between anisotropic monolayer graphene and bi- 
layer graphene in the presence of a translation between 
the layers. This report is organized as follows: In the next 
section, we introduce the tight-binding Hamiltonian for 
a monolayer and consider its limits when ts > t/3 and 
ts < t/2 corresponding respectively to the emergence of 
new Dirac points in the vicinity of the three inequivalent 
M points of the reciprocal space and to their merging 
at the K points. In section III, we recall the effective 
low energy Hamiltonian for bilayer graphene and we dis- 
cuss its relation with the monolayer Hamiltonian, with 
a correspondance between the relevant parameters. We 
briefly conclude in the last section and we stress that a 
simple model with a single parameter ts/t allows for the 
continuous description of the evolution between different 
scenarios of merging (or emergence) of Dirac points with 
opposite or same Berry phase. These remarks provide a 
simple tool to describe in a coherent picture a large range 
of different physical situations. 



Let us first recall the tight-binding Hamiltonian for the 
monolayer honeycomb lattice : 

W-rnono = -t^2(albj + h.c.)-t 3 ^ (a\bj + h.c.) . (1) 

(iJ) «<M»> 

where ai (aj) annihilates (creates) and electron on a A 
site and there is a similar definition of the bi operators 
on the B sites (Fig. 1). The first term describes the 
coupling between nearest neighbors (t being the hopping 
amplitude), and the second term describes the coupling 
between third nearest neighbors (ts being the hopping 
amplitude). These two terms correspond to coupling be- 
tween sites of different sublattices. We have not consid- 
ered here the coupling between second nearest neighbors 
(which couples sites of the same sublattice), because it 
is not relevant here for our purpose. This Hamiltonian 
can be rewritten in the form H = ^ with 
Wc = («k4) and 4 



H-monoO^) 



t/(k) + t 3 / 3 (k) 

t/*(k) + i 3 / 3 *(k) 



where the functions /(k) and /^(k) are given by 



/(k) 



(2) 



(3) 



/ 3 (k) = e - 2ik -*i + e -2*k.<5 2 + e -2ik.S 3 (4 ) 

and the vectors Si connect one atom to its three nearest 
neighbors (Fig. 1): 

s i= a (^Y>\)> <*2 = a(-^,^), <5 3 = a(0,-l), 

(5) 

a being the interatomic distance. 

In graphene, the ts term is small. 5 However, it is of in- 
terest to imagine a larger value of this parameter because 
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FIG. 1. Honeycomb lattice. The Si connect one atom to its 
nearest neighbors. The are Bravais lattice vectors. 



it has a quite interesting effect on the evolution of the 
spectrum, as has been theoretically considered in ref.[2]. 
When ts increases and reaches the critical value t/S, a 
new pair of Dirac points emerges from each of the three 
inequivalent M-point of the reciprocal space (see Fig. 
2.a,b), following the universal scenario predicted for the 
apparition of Dirac points in a 2D crystal in the vicinity 
of time-reversal symmetry points. 6,7 Writing k = \lm +q 
, we obtain the following " universal Hamiltonian" in the 
vicinity of £3 = t/S (keeping the leading order terms) 



H-univ. (q) 
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where q = (q x ,q y ) and where the parameters m*,c, A* 
can be related to the original band parameters. Here we 
find A* = t — 3^3, c = 2t and m* = 2/t. The parameter 
A*, when becoming negative (ts > t/S), drives the emer- 
gence of a new pair of Dirac points at the M-point (Fig. 
2.b). The distance between the new Dirac points is given 
by Aq 2V-2m*A* = Ay/St$/t - 1 

When increasing further £3, the new Dirac points ap- 
proach the K and K' points, so that each initial Dirac 
point sitting at the ) points is surrounded by three 
Dirac points (with opposite Berry phases, see Fig. 2. a). 8 
These Dirac points merge at the critical value ts = t/2, 
and the spectrum becomes quadratic around K^, as al- 
ready noticed in Ref. 2 . We will show in the next section 
that there is an exact correspondance between the Hamil- 
tonian (1) with ts =t/2 and the low energy Hamiltonian 
for bilayer graphene. When ts — t/2, 1-Lmono takes a 
new universal form (keeping leading order terms) in the 
vicinity of the K' point 



V 



2m* 



cq ] 



2m* 



cq + A 







(7) 

where q — q x + iq y (The Hamiltonian in the vicinity of 
the K point is obtained by the substitution q — >• —q^). 






FIG. 2. Iso-energy lines in the vicinity of the ^ and M 
points, for different values of the parameter £3. (a) ts — 0, 
(b) £ 3 = 0.35t, (c) t 3 = 0.40£, (d) t 3 = 0.5t, , (e) t 3 = 0.65£. 
The vicinity of the Dirac points is indicated in red, as well as 
their associated Berry phase (within a factor tt). 



As we shall see in the next section, this is precisely 
the low energy effective Hamiltonian for sliding bilayer 
graphene. 9 12 Starting from (2), we find m* = 4/9£ and 
c = 3(ts — t/2) . Here the parameter A = 0. We comment 
on its effect in the next section. We see that the param- 
eter c controls the number of Dirac points. When c = 0, 
the spectrum is quadratic. When c becomes finite, three 
Dirac points emerge in a trigonal arrangement. 3 The ori- 
entation of this arrangement depends on the sign of m*c, 
that is the sign of t(2ts — t). 




FIG. 3. Iso-energy lines for the Hamiltonian (7), here with 
A = 0. When c = the spectrum is quadratic. When c / 0, 
the spectrum splits into four Dirac points, keeping conserved 
the total Berry phase. The orientation of the Dirac points 
depends on the sign of m*c. In changes when £3 becomes 
larger than t/2. These figures (here for m* > 0) have to be 
compared with the spectrum in the vicinity of the K/ point 
shown in Figs. 2.c-e for £3 < t/2, £3 = t/2 and t$ > t/2. 
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III. BILAYER 

We now consider the Hamiltonian describing a bilayer 
of honeycomb lattices arranged in the so-called Bernal 
(or A — B) stacking (Fig. 4. This Hamiltonian reads 1 



Here /k is a short notation for /(k). The correspond- 
ing energy spectrum exhibits four bands. At low energy, 
two bands are located around energies ±71, while the 
two other bands touch each other at zero energy. When 
7o|/k| <^ 7ij the low energy Hamiltonian describing the 
two lowest bands reduces to a 2 x 2 matrix 



Km = - 7o ^2(4 b 3 + h - c -) 

(id) 

-7o$Z(a|&j + h.c.) 

(id) 

~ 7i ^2(a\bj + h.c.) 

3 

-73^(4^ + h.c.) 



(8) 
(9) 



The a,i,bi operate on one layer (A and B sites) and the 
(ii,bi operate on the other layer. The 70 terms describe 
the coupling between nearest neighbors in each layer. 
The 71 term couples sites from different layers which are 
on top of each other, the A and B sites in Fig. 4. The 73 
term describes the coupling between A and B sites be- 
longing to different layers. In bilayer graphene, one has 
approximatively 73 c± 71 ~ 70/ 10. 




FIG. 4. Graphene bilayer in the Bernal stacking. The atoms 
A and B are on top of each other. The pair (A-B) is denoted 
C. This picture shows that the indirect coupling between A 
and B sites results from second order processes where the 
electron jumps from A to C — (A, B) and then from C to 5. 
There are two different processes indicated on the figure, one 
with black full lines, one with a dotted line. The first process 
implies to possible paths and its amplitude is thus twice larger 
than the second one. 



nil 
rth 







2 

70 r*2 

\ 7i 



— /k +73/k 

71 



73/k 







(11) 



J 



and the electron resides mainly on the sites A and B 
which do not face each other (denoted C on Fig. 4). As- 
suming 73 = 0, the spectrum is quadratic in the vicinity 
of the K( ) points. Due to a finite 73 the finite structure 
of the low energy spectrum (trigonal warping) exhibits 
four Dirac points with a linear dispersion. 3 

We now discuss the relation between the Hamiltonians 
(2) and (11). In order to have a simple intuitive picture, 
it is instructive to discuss the expansion of the function 



2ik.<5! 



3 2ik.<5 2 



3 2ik.(5 3 



_|_ 2 ^ e *k.(<5i+<5 2 ) _|_ e ik. (<5 2 +<5 3 ) _|_ e «k. (61+63)^ ^2) 

This structure can be interpreted physically as follows : 

2 ~ 

The effective hopping term between A and B sites 
results from an indirect second order hopping via C = 
A, B) sites : an electron hops from A to C with an en- 
ergy cost 71 on site C and then hops from C to B). Fig. 
4 shows the two different possible types of paths, one via 
two intermediate steps at the corners of a losange (asso- 
ciated to a translation vector —£2 — £3 in Fig. 4), the 
other one along a straight line, associated with the trans- 
lation vector —2^i in Fig. 4). Since there are twice more 
processes of the first kind, we immediately understand 
the structure (12) of the function f£. 

Moreover, since Si + £2 + £3 = 0, the function f£. can 
be rewritten simply as 



/(k) 2 = / 3 *(k) + 2/*(k) 



(13) 



so that we get the simple correspondance between 
the monolayer Hamiltonian and the low energy bilayer 
Hamiltonian : 



This Hamiltonian can be rewritten in the form T~L = 



Ek4^( k )^ with 
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(10) 



% l ™(k) — H^ ono (k) 
with the correspondance 



t = 73 + 2 — , £3 = — 
7i 7i 



(14) 



(15) 
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FIG. 5. Parameters m* and c of the universal Hamiltonian 
^i'univ (V an ^ their relation to the tight binding models for 
monolayer and bilayer lattices. E s — m*c 2 /2 is the energy of 
the saddle point between the three satellite Dirac points and 
the central one. 



that is 73 = t — 2ts. 

In the limit where 73 is small, or alternatively when 
ts — t/2, both Hamiltonian reduce to the form l~H univ 
written in (7), with the parameters given in the table 
5. It describes a set of three Dirac points arranged at 
a distance Aq = 2m* c of the central one and separated 
by a saddle point (Van Hove singularity) at energy E s = 
m*c 2 /2. 

Finally, there has been several works on deformed bi- 
layer graphene, described with the Hamiltonian H l £ w — 



AI, where the parameter A is finite and is related to 
a deformation (translation or rotation) between the lay- 
ers. This is precisely the Hamiltonian (1) written for the 
monolayer. One could wonder what is the significance 
of the parameter A for the monolayer Hamiltonian. The 
answer is simple : it is obtained by a uniaxial anisotropy 
of the hopping integrals between first nearest neighbors. 
Assuming that the hopping integral is different t' 7^ t 
(and also t' 3 7^3) along the vertical axis, we obtain that 
the parameter A is given by A = t — t' + £3 — t 3 . 



IV. CONCLUSION 

We have shown a simple correspondance between the 
Hamiltonian of monolayer with uniaxial anisotropy and 
a large nearest neighbors coupling £3, and the low en- 
ergy Hamiltonian for distorted bilayer graphene. This 
provides a simple tool to study apparently different sys- 
tems. For example, the monolayer Hamiltonian could be 
used to study the edge states of the bilayer graphene. 
In addition it is noticeable that a single tight-binding 
Hamiltonian with two parameters ts/t and t' jt describes 
continuously the evolution between different scenarios of 
merging or apparition of Dirac points. 
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